Hopping charge transport in organic materials 
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General properties of the transport of charge carriers (electrons and holes) in disordered organic 
materials are discussed. It was demonstrated that the dominant part of the total energetic disorder 
in organic material is usually provided by the electrostatic disorder, generated by randomly located 
and oriented dipoles and quadrupoles. For this reason this disorder is strongly spatially correlated. 
Spatial correlation directly governs the field dependence of the carrier drift mobility. Shape of the 
current transients, which is of primary importance for a correct determination of the carrier mobility, 
is considered. A notable feature of the electrostatic disorder is its modification in the vicinity of the 
electrode, and this modification takes place without modification of the structure of the material. 
It is shown how this phenomenon affects characteristics of the charge injection. We consider also 
effect of inter-charge interaction on charge transport. 
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I. INTRODUCTION 

Charge carrier transport in organic materials is a vi- 
brant area of research for over three decades. To a large 
degree this attention is motivated by the constantly grow- 
ing application of organic materials in various electronic 
devices. In this paper we will completely put aside this 
aspect of the organic electronics and focus entirely on 
the fundamental properties of the conductivity of organic 
materials. Reviews of advances in the application of elec- 
tronic organic devices may be found in recent books . 
In addition, we do not consider properties of highly con- 
ductive organic polymers (doped polyaniline, polyacety- 
lene and others) and limit our attention to more tra- 
ditional materials having low conductivity. By the usual 
classification they may be considered as amorphous semi- 
conductors with wide bandgap. These materials are poly- 
meric or low molecular weight organic glasses. Typical 
example of the conducting organic glass is a material cre- 
ated by doping of the inert polymer binder (polycarbon- 
ate or polystyrene) with molecules, providing conducting 
properties (mostly, aromatic amines, hydrozones, or ni- 
triles); mass fraction of the dopant in most cases falls in 
the range 30%-100% Q. 

Most important features of the charge transport di- 
rectly follow from the basic structural features of or- 
ganic glasses. These glasses are molecular materials with 
rather weak interactions between molecules and, at the 
same time, they have significant disorder in positions and 
orientations of molecules. This means that all relevant 
states are localized and charge carrier transport occurs 
by the hopping mechanism 5]. It is well known that 
for the hopping transport the disorder in the material 
is of paramount importance. Traditionally, total disor- 
der in amorphous materials is subdivided into energetic 
and positional disorders. The energetic disorder includes 
random fluctuations of the positions of energetic levels 
of molecules, while positional disorder describes fluctu- 
ations of the positions of molecules that affect hopping 



rates between transport sites but do not lead to the fluc- 
tuations of energies. There is a general agreement in the 
community that the energetic disorder is of primary im- 
portance for the hopping transport (at least for the case 
where the concentration of transport sites is not too low) 
6]. In this paper we consider general statistical proper- 
ties of the energetic disorder in organic glasses, how they 
affect various features of the hopping transport in organic 
materials, and provide comparison with well-established 
experimental facts. 



II. BASICS OF HOPPING TRANSPORT IN 
AMORPHOUS MATERIALS: DENSITY OF 
STATES 

Typically, a hopping model is constructed as a reg- 
ular cubic lattice with the lattice scale a, its sites are 
randomly occupied by transport molecules with the oc- 
cupation fraction c, and a random energy Ui is assigned 
to each site. In most cases the hopping rate has a Miller- 
Abraham (MA) form 'f\ 



Pi-^.j = t^oCxp(-27ry) 




* lnovikov@elchem.ac.rul 



(1) 

where is some characteristic frequency of hops, = 
\rj — Til, and 7 is a wave function decay parameter for 
transport sites; in organic materials 7a ~ 5 — 10 If 
an electric field E is applied, then the random energy in- 
cludes an additional term —efiE. Distribution of random 
energies is typically considered having a Gaussian Q 
or exponential form (loj . In this paper we limit our con- 
sideration to the Gaussian density of states (DOS). This 
particular form is widely recognized as more suitable for 
description of organic amorphous solids ^ . Apart from 
the experimental evidence, the Gaussian DOS naturally 
arises in simple but reasonable models of disordered or- 
ganic materials. For example, the simplest model of the 
polar glass is the lattice with sites occupied by randomly 
oriented dipoles [HI, [lH . If fraction of occupied sites is 
not too low (c ~ 1), then the resulting DOS P{U) has a 
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Gaussian form with the width 



2.35 



epc 



,1/2 
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here p is the dipole moment, and e is the dielectric con- 
stant of the medium [l2| (numeric coefhcient 2.35 is spe- 
cific for the simple cubic lattice). For typical values of p, 
£, and c the magnitude of disorder a ~ 0.05 — 0.1 eV. For 
c <C 1 there is an intermediate asymptotics 



P{U) cx 
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but in practice this regime cannot be observed due to 
inevitable additional contribution to the DOS from other 
sources of disorder [T2| . 

Particular shape of DOS directly governs the mobility 
temperature dependence /i(T) (for a low field region). 
Gaussian DOS leads to 



In jj, '. 



\kT) 
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where A ~ 1 is some constant, and this particular de- 
pendence is believed to be the most pro perly suited for 
description of experimental data [13l - [l5| . Roots of Eq. 
([4]) may be understood if we calculate the position of the 
maximum of the density of occupied states 



Pocc{U) oc P{U) exp(-t//fcT) cx exp 



2a^ 



U 
kT 
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in the Gaussian DOS, and then use ?7max as an equivalent 
of the temperature-dependent activation energy. 

Motivated by this evidence, Bassler suggested a Gaus- 
sian Disorder Model (GDM) for the description of hop- 
ping charge transport in amorphous organic materials [6| . 
This model was claimed to be a universal model suit- 
able to describe charge transport in any disordered or- 
ganic material. Main ingredients of the GDM are the as- 
sumption of the Gaussian DOS and validity of the Miller- 
Abrahams hopping rate ([T]) . It was assumed also that the 
distribution of random energies has no spatial correlation 
at all; the correlation function is {U{ri)U{rj)) — a^Sij. 
The major part of results has been obtained for the GDM 
by means of Monte Carlo simulation. For example, it was 
found that for the GDM A w 4/9 in Eq. 

Another important parameter affecting mobility is an 
applied electric field E. For a long time it was known that 
organic glasses do demonstrate a strong field dependence 
of the mobility [IB] , and the so called Poole- Frenkel 
(PF) dependence 



In /i cx const + \^ 



(6) 



is usually a good description of the experimental data. 
While the GDM fairy well describes the mobility tem- 
perature dependence, it totally fails in describing its field 
dependence. There are claims that the GDM do repro- 
duce the PF dependence Eq. ©, though in a limited 
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FIG. 1. Field dependent mobility in the GDM for differ- 
ent values of kT/o (from the top curve downward); /io = 
ea'^i/o/cr exp (— 27a). Straight lines show best fits for the lin- 
ear regions of the curves. 



field range |6|. This statement is not true. Field depen- 
dence of /i in GDM for low and moderate fields is much 
better described by a linear dependence 



1,17(^-2,05 
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This equation summarize the results of the Monte Carlo 
simulation (see Fig. [1]), more thoroughly described else- 
where [l3l. In fact, the parabolicity of the plot In// vs 
i?^/^ is clearly seen even in Fig. 7 of the Bassler's paper 
0. Close inspection of Fig. [T] (as well as Fig. 7 of pa- 
per ^) shows that the only field range that can mimic 
the PF dependence (O is the region eaE/a > 1, where 
the mobility curve, plotted as In /x - E, begins to deviate 
from the straight line. The reason for this deviation is 
the use of the particular hopping rate, that is the MA 
rate: it has a property that for very strong electric field, 
where almost all hops occur downward in energy, carrier 
velocity is saturated and, hence, the mobility begins to 
decay with E as 1/ E. By its very nature, this mechanism 
cannot provide a good linearity of the dependence In vs 
£;i/2 fQj. ^}ig range 10^ — 10® V/cm, where it is routinely 
observed (if cr = 0.1 eV, then for the typical scale 
a = 1 nm eaE/a w 1 for = 1 x 10® V/cm). In some 
materials the PF dependence was tested even for much 
wider range, such as 8 x 10^ — 2 x 10® V/cm [l8|. 

Inability to describe the mobility field dependence in- 
dicates that some important element is missing in the 
GDM. We are going to find this missing element. 
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FIG. 2. Distribution of site energies U in the lattice model of 
dipolar glass. A sample with the size of 50 x 50 x 50 lattice 
sites is shown. Black and white spheres represent the sites 
with positive and negative values of U, correspondingly, while 
the radius of a sphere is proportional to the absolute value of 
U. Sites with small absolute values of \U\ (less than a) are 
not shown for the sake of clarity. 



III. ORGANIC GLASSES: KINGDOM OF 
SPATIAL CORRELATIONS 



Long range spatial correlations of the energy 
landscape in organic glasses 



The GDM is based on the assumption of the non- 
correlated nature of the random energy landscape in 
organic materials. It turned out that this assumption 
proved to be spectacularly wrong in organic glasses. In- 
deed, we already mentioned that organic glasses have 
very low concentration of free intrinsic carriers and, 
hence, do not provide effective screening of electrostatic 
interactions. At the same time, they have high concen- 
tration of permanent dipoles and quadrupoles. These 
molecules provide long range (unscreened) random con- 
tributions to the total random energy of charge carri- 
ers. Sum of random terms, slowly decaying in space, 
inevitably produce strongly correlated random energy 
landscape. For dipoles the corresponding correlation 
function decays as 



(8) 



while for quadrupoles [2C 

C,(f) p 



(9) 



Dimensionless parameters Ad and Aq equal to 0.76 and 
0.5, correspondingly, for the simple cubic lattice [13, [HI- 
Quite probably, the correlation functions of dipolar or 
quadrupolar type are the most common ones in organic 
glasses, because the model of randomly oriented and 
located dipoles provides a good approximation for the 
polar organic glasses, and the corresponding model of 
quadrupoles gives a reasonable model for non-polar or- 
ganic glasses. 

Long range correlations mean that the random energy 
landscape in organic glasses has a natural cluster struc- 
ture: sites with close values of U tend to group together 
(see Fig. [2]). Abundance of clusters in organic glasses 
may be characterized by the cluster distribution on size: 
for the dipolar glass the asymptotics for the number of 
clusters Ug having s sites decays as 



Inn, 



oc 



a/3 



(10) 



while for the non-correlated GDM 



lung c>c — s, (11) 
which gives a huge difference for s ^ 1 [2^ [2^ . 



B. How correlations affect charge mobility 
dependence 

Correlation nature of the energy landscape C/(r) di- 
rectly dictates major features of the charge transport and 
injection in organic glasses. For example, the field depen- 
dence of the quasi-equilibrium mobility can be under- 
stood from the following simple consideration (for more 
thorough consideration see Ref. [l^l). Suppose that the 
carrier is located at the bottom of the potential well with 
the energy U{0). Mobility is determined by the typical 
time for the carrier to reach a saddle point with the en- 
ergy U{r) — eEf, where r is the distance to the saddle 
point from the bottom of the well. That time can be 
estimated as 



t ~ to exp 



U{r) - U{0) - eEf 
kf 



(12) 



where to — 1/t'o, and the average time for the Gaussian 
random landscape U{r) is 



(i) ~ io ( cxp 



U{f} - U{0) 



kT 



exp 



to exp ■ 
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'kT 



to exp 
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C(0) - C{f) eEf 
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Assuming charge transport in the dipolar glass (DG) with 
the correlation function , we can calculate the critical 
size of the potential well that provides the maximal es- 
cape time 



d{t) 



-0, 



aAd 



dr ' V eEkT 

and the mobility is estimated as 



1/2 



(14) 



^ cull {t (rcr)) oc exp 



f-V 

\kT) 



kf \ kT 



(15) 

This result provides a leading asymptotics of the exact 
solution of ID transport problem ,25,^]. Typical mag- 
nitude of the slope of the mobility field dependence In/i 
vs E^/^, calculated using Eq. ([T5|). agrees well with the 
estimations that follow from the experimental data for 
a, obtained from the low field mobility temperature de- 
pendence [3]. This means that in polar materials the 
dipolar a (estimated from the mobility field dependence) 
provides a dominant part of the total a (estimated from 
the mobility temperature dependence). Extensive com- 
puter simulation of the 3D transport generally confirms 
Eq. (lisp and only modifies numeric parameters in this 
relation. Results of the simulation may be summarizes 
as a phenomenological relation 



lnii = - 
Mo 











KkTJ 



(16) 

where Co ~ 0.78, and T w 2 [21]. 

In a more general case of the algebraic correlation func- 
tion 



C(r) = A„a2(^)' 



(17) 



the mobility in ID case and strong disorder a/kT ^ 1 is 
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For the GDM formally n — >■ oo and the result agrees well 
with the leading field-dependent term in Eq. Q. 

This analysis indicates not only that the GDM fails to 
capture a very important characteristic of organic glasses 
(spatial correlations), but also that there is no unified 
mobility field dependence for different classes of organic 
materials. The functional form of the mobility field de- 
pendence depend on the spatial decay of the correlation 
function and is, indeed, different in polar and non-polar 
organic glasses [1^. This conclusion of ID analysis is in 
good agreement with the result of 3D simulation for the 
model of quadrupolar glass (QG), serving as a reasonable 
model of the non-polar organic glass (Fig. [3]). Simulation 



FIG. 3. Mobility field dependence in the QG model for differ- 
ent values of kT/o (from the top curve downward); straight 
lines in (a) indicate the fit for Eq. (|19p . Plot of the sim- 
ulation data in the usual PF presentation In /i vs E^^'^ (b) 
demonstrates much stronger deviation from the linearity in 
the weak field region (straight lines serve as a guide for an 
eye). 



data can be reasonably well described by the equation 



Q 



with Cq « 0.87 and Tq « 1.91 [Sfl]. 

For the non-correlated Gaussian landscape mobility is 
controlled by the carrier release from the deep states to 
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the neighbor sites having higher energy, hence the shift 
of the carrier energy in the apphed field lead to the linear 
field dependence 



In /i cx eaE/kT, 



(20) 



which again agrees well with Eq. ([7]) and the limiting 
case of Eq. (fTSj). 

An interesting difference between correlated energy 
landscapes and the GDM is a non-existence of the so- 
called transport energy in organic glasses [3l| . This con- 
ception was extensively used for an analysis of the hop- 
ping transport in amorphous materials 32- 34]. Trans- 
port energy serves as an analogue of the mobility thresh- 
old for the hopping transport. It emerges as a result of 
the competition of two opposite tendencies for a carrier, 
escaping from low energy sites: fast decay of the hopping 
rate Eq. ([T]) with distance facilitates hops to the nearest 
neighbors, while the the probability to find a neighbor 
with not too high energy (and, hence, not too low hop- 
ping probability) increases with distance. As a result, 
in the GDM there is an optimal hopping distance and 
the optimal final energy (transport energy) for a carrier, 
which does not depend on the initial carrier energy. In 
the correlated landscape site, close in space, are the sites, 
close in energy, too, and the optimal final energy does 
not exist (the most probable final energy shifts with the 
variation of the initial carrier energy [31|, see Fig. H]). 




FIG. 4. Position of the maximum of distribution of final car- 
rier energy Uf in DG (■), QG (a), and GDM (□), corre- 
spondingly, in dependence on the initial carrier energy Ui for 
kT/a = 0.25. Lines are provided for the convenience. 

The most important conclusion of this section is that 
the functional form of the mobility field dependence is 
directly dictated by spatial correlations in the particular 



organic material. The hope to provide a single universal 
formula for fJ.(T, E) is futile, and the major watershed 
divides the polar and non-polar organic materials. At the 
same time, analysis of experimental data indicates that 
a reliable determination of the mobility field dependence 
is not a trivial task. For all non-polar organic glasses 
rather limited experimentally tested field range is typical 
(one order of magnitude or even more narrow i35l - 
|38|). Taking into account inevitable experimental errors, 
this means that the reliable discrimination between the 
PF dependence ^ and the true quadrupolar dependence 
is hardly possible (a good example is discussed in the 
paper [2_g]). 



IV. 



CURRENT TRANSIENTS: CONCEPTION 
OF UNIVERSALITY 



A most important method for the experimental deter- 
mination of the hopping mobility in organic materials is 
a time-of- flight (TOF) method |4|. Here the experimen- 
tal sample is a slab of organic material with thickness L 
sandwiched between two electrodes. Typically, electron- 
hole pairs are generated in the vicinity of one electrode 
(though a variant with the uniform bulk generation is 
used, too [3^, but for this variant some serious difficul- 
ties in the interpretation of the experimental data take 
place [40]), carriers of one sign are instantly absorbed 
by the electrode, and carriers of the opposite sign drift 
through the slab to the collecting electrode under the ac- 
tion of the applied field E. In the TOF method the raw 
experimental data is the time dependence of the current 
/(t), generated by drifting carriers (current transient). 

In most cases transient demonstrates an initial short 
spike, indicating a spatial and energetic relaxation of car- 
riers, then a plateau (at this stage carriers move with the 
almost constant velocity v), and finally a fast decay, in- 
dicating an arrival of carriers to the collecting electrode. 
These are typical features of so-called non-dispersive (or 
quasi-equilibrium) transport, and the most important 
property is the independence of a mean velocity (v) on 
L, well established in experiments [4lj. Mobility is cal- 
culated as 



M = 



E 



L 



(21) 



where td is some characteristic drift time. In most pa- 
pers td was chosen as a time tg of the intersection of the 
asymptotes to the plateau and tail of the current, while 
in some papers td is the time of the current to decay 
to one half of the plateau value. 

In some cases (mostly for low temperature) plateau of 
the transient is not well defined and current demonstrates 
monotonous decay. This is the case of the dispersive 
(non-equilibrium) transport and the shape of the tran- 
sient is analyzed in double logarithmic coordinates In / 
vs Int according to the popular Scher-MontroU model 
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m 



t<tT, 
t>tT, 



(22) 



here tr oc L^/" is some characteristic time and < a < 1 
is the dispersive parameter. In this regime the mean 
velocity depends on L as (v) oc L/tx oc L^~^/". 

If we discuss charge transport in term of mobihty, then 
all we need from the experiment is the drift time td (be 
it tQ, ti^2i or tx)- Naturally, the whole shape of the 
transient could provide additional valuable information 
about the transport mechanism. In fact, sometimes even 
a proper discussion of the mobility field dependence re- 
quires a clear distinction of the mobilities, obtained with 
to and ti/2, i.e. an explicit use of some information about 
shape of the transient [431] . This remark explains the sig- 
nificance of the parameter 



W 



tl/2 — to 



1/2 



(23) 



which provides a simplest robust integral characteristic of 
the shape. Other important features of the transients are 
short and (especially) long time asymptotics (the initial 
spike and tail of the transient). 

All theoretical approaches indicate that for the Gaus- 
sian DOS a constant (independent of time) velocity 
V = (v) does eventually emerge and the quasi-equilibrium 
steady state is achieved [l]; for this reason in the pa- 
per we limit our consideration to quasi-equilibrium non- 
dispersive transients. Experimental studies show that a 
general picture could be summarized in the following way. 
For high temperature a/kT < 3 parameter W decreases 
with L as oc [11, Ell, which is a fingerprint 

of the classical diffusion. Indeed, it was found that the 
usual classical diffusion is a good approximation for a 
description of the shape of transients in this situation 
[4^147.]. 

For lower temperatures current transients have a dif- 
ferent form [H, |4^: for such transients W{L) fa const. 
Hence, if the time scale is properly re-scaled, then tran- 
sients for different L approximately collapse to a single 
universal curve (for this reason this phenomenon was 
dubbed universality) . It was found that universality with 
respect to L usually means the universality with respect 
to E, but not T [481]. 

Previously, the very conception of universality has 
been invariably attributed with the dispersive transport. 
Schein et al. (48^] specifically emphasized that their's data 
suggest a universal behavior for the non-dispersive trans- 
port with the mean carrier velocity, which does not de- 
pend on L. This fact clearly indicates that a common 
perception (see, e.g. j46|) that a transport with the well- 
defined constant velocity is invariably the transport, de- 
scribed by the usual diffusion equation, is certainly not 
true for organic glasses. In fact, this phenomenon has 
been well established long ago for some models (see the 



excellent review [49|), but did not attract much attention 
in the transport community. 

Unfortunately, we have no reliable analytic results for 
the shape of transients in the random medium in 3D case. 
Most interesting exact results have been obtained for ID 
transport [sO] ■ This paper was specifically devoted to the 
study of the diffusing coefficient D in correlated Gaussian 
energy landscape with DG correlation function Eq. ([S|). 
It was found that the Einstein relation 

likT 



D = 



does not hold, but the modified Einstein relation 



D 



kT d (v) 
~e dE 



(24) 



(25) 



is valid, which transforms to Eq. (f24t for /i(-E) = const. 
Yet the major result of the paper [5Q] is the very exis- 
tence of the diffusion coefficient for all T and i? in ID 
transport. This is a drastic contradiction with the uni- 
versality, found in experimental papers 45, 4S] ■ 

Monte Carlo simulation agrees well with the experi- 
ments. Indeed, for high temperature the transients are 
diffusive [5l| . while for the low temperature a univer- 
sal behavior emerges (Fig. [5]). Note, that in our case 
the transport is definitely non-dispersive because a single 
curve is produced by re-scaling the time as t ^ t (v) /L 
and velocity as v ^ v/ (v) (the later one is needed only 
for checking of the i?-universality). Is the universality 
a distinct feature of the 3D transport or true diffusive 
behavior still emerges in 3D case for a very long time, re- 
mains an open question. It is worth to note that accord- 
ing to the simulation data, the modified Einstein relation 
Eq. (|25p does not hold for the 3D transport in the DG 
model even for high temperature [5l| . 



V. ENERGETIC DISORDER NEAR THE 
ELECTRODE 

We already noted that the total energetic disorder in 
organic glasses is mainly the electrostatic disorder, pro- 
vide by randomly located dipoles and quadrupoles. For 
this reason, it demonstrates a very unusual property: 
modification of the statistical properties of the disorder 
near the conducting electrode without any modification 
of the structure of the material. 

In fact, in any case structure of the organic material in 
the vicinity of the electrode should be quite different from 
the bulk structure. We should expect different packing 
of spacious organic molecules, accumulation of impuri- 
ties, partial degradation of the organic material etc. For 
all these reasons the energetic disorder at the electrode 
is very different from the bulk disorder (typically, it is 
greater than the bulk disorder). 

Yet in organic glasses there is an opposite general con- 
tribution, leading to the decrease of the electrostatic dis- 
order at the electrode. The electrostatic energetic disor- 
der is directly proportional U{r) = eip{r) to the disorder 
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FIG. 5. Test of the universality of current transients for the 
DG model; kT/a = 0.19 and c = 0.3. Carrier velocity v is 
directly proportional to I{t). a) Transients for the electric 
field in the range 0.02 - 1.2 eaE/a and L = 20,000a. b) 
Transients for eaE/a = 0.02 and L varying from 2,000 to 
20,000 lattice scales. 



in the spatial distribution of electrostatic potential f{r), 
generated by molecular dipoles or quadrupoles. In or- 
ganic layers sandwiched between conducting electrodes 
the potential ip{r) must obey a boundary condition at 
the electrode surface: here potential should be a con- 
stant. Thus, at the electrode there is no electrostatic 
disorder at all, irrespectively to how disordered is the 
material in the bulk. This means that the magnitude of 
the dipolar or quadrupolar disorder increases while going 
away from the electrode, asymptotically reaching its bulk 
value. 

Some decrease of the disorder at the surface of organic 
material is inevitable for any model of organic glass (just 
because there are more neighbor molecules in the bulk of 
the material) , but the magnitude of the effect in the case 
of electrostatic disorder is much greater than in the case 
of short range interactions. For example, for the sim- 
ple model of the interaction with the nearest neighbors 
only and simple cubic lattice we have cr^^j.^^^^ 
while for the dipolar disorder CTg^^rface 
disorder depends on the distance z from the electrode as 



0.3cr^^jj^ and rms 



2 

"bulk 



1 



Oo 

2z 



ao = Ada 



(26) 

(here (Tsurfacc is the magnitude of the disorder in the first 
layer of organic material, directly adjacent to the elec- 
trode). Spatial correlations at the electrode do differ too; 
this phenomenon has no analogue for the short range dis- 
order. A direct calculation of the correlation function 
C(r) near the electrode gives 



1 



1 



(27) 



where rj^ — p"^ + {zi ± 2:2)^ and p is a 2D vector oriented 
along the electrode plane [H^ . Hence, at the electrode the 
dipolar glass is much less correlated in comparison with 
the bulk (Fig. C{zi,Z2,p) oc ziZ2/p^ for p > 21,^2, 
and clusters are elongated perpendicular to the electrode 
plane. 

Decrease of the disorder at the electrode and change 
of the spatial behavior of the correlation function have 
very important implications for the charge injection. In 
the absence of these phenomena, injection current in or- 
ganic glasses demonstrates formation of channels where 
the current density is much greater than the average den- 
sity Such channels originate from particular spots 
at the electrode, where clusters of sites with low energy 
facilitates injection. Reduction of the disorder at the 
electrode and modification of the spatial behavior of the 
correlation function lead to the more uniform distribu- 
tion of the injection current over the electrode and dra- 
matically reduce current channeling. This decreases local 
overheating in a device and improves its performance. 

In conclusion, we would like to emphasize that by no 
means the presented calculation gives a full and accurate 
estimation of the drop of a at the electrode. In fact, as it 
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FIG. 6. Distribution of the carrier random energy (7(r) in the 
bulk of the organic material (a) and in the organic layer closest 
to the electrode (b). All symbols are exactly the same, as used 
in Fig[2] Note the significant decrease of the amplitude of the 
disorder at the electrode and its much less correlated nature. 
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FIG. 7. Dependence of the effective temperature on the car- 
rier density p for kT/a = 0.3 (■) and kT/a = 0.4 (□) for 
the DG model (a) and GDM (b), correspondingly. Insets 
show the curves ij.{E) for kT/a — 0.3 and different p: 0.0016, 
0.008, 0.016, and 0.032, from the bottom curve upward, cor- 
respondingly. For the DG T^s was calculated by fitting the 
mobility data to Eq. (|16|) . while for the GDM the mobility 
data was fitted to Eq. ((T)). Lines are provided as guides for 
an eye. 



was mentioned already, there are many other reasons for 
the variation of magnitude of energetic disorder at the 
electrode and in many cases the total disorder increases 
near the electrode. There are three major reasons to 
consider the particular drop of the electrostatic disorder 
at the electrode. First, this effect provides a significant 



contribution to the variation of the magnitude of total 
disorder, in polar materials it could reach about 50% of 
Cbuik- Second, this effect is universal, it does not depend 
on particular details of the structure of the material at 
the electrode (though, of course, an actual magnitude of 
the drop does depend on the structure of the organic glass 
near the electrode). And third, it is very unusual, we al- 
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ready noted that for many other contributions we should 
expect the increase of the disorder at the electrode. 

VI. TRANSPORT OF INTERACTING 
CARRIERS 

A. Dynamic effects of interaction 

If charge carrier density is not very low, we cannot ne- 
glect the Coulomb interaction between carriers. The rel- 
evant interaction strength parameter is Uc — e^n}^^ / ea, 
where e ~ 2.5 — 3 in organic materials and n is a carrier 
concentration. The maximal value of Uc for n ~ l/o? 
in organic glasses is U^^^ ~ 5. This means that for 
high concentration of carriers we cannot neglect the inter- 
charge interaction. 

Theoretical studies of the effect of carrier density are 
not numerous [ssl . ISGj . High density of carriers could 
affect drift mobility in opposite ways. Small fraction of 
carriers could occupy deep states, thus providing a pos- 
sibility for remaining carriers to avoid trapping and ac- 
quire much higher mobility. At the same time, charge- 
charge interactions could provide an additional energetic 
disorder in the material. This is indeed the case for the 
simplest model where all charges except one are immov- 
able [13, HI]. In this case, the greater is the density of 
static charges (i.e. energetic disorder), the smaller is the 
mobility. 

Usually, in theoretical studies the mean field approxi- 
mation has been used and charge-charge interaction has 
been totally neglected [1, HI]. This means that these 
studies dealt only with the effect of filling of deep states: 
it is assumed that the effects of interaction could be later 
effectively included via the mobility dependence on the 

mean local electric field (^E\o^, which in turn is con- 
nected to the mean local charge density p by the Poisson 
equation 

div(4oc)-%. (28) 

This line of reasoning totally neglects dynamic cor- 
relations. In addition, quite frequently true quasi- 
equilibrium mobility is formed by the averaging over large 
domains of the disordered material (see, for example, Ref. 
0). In this situation the very conception of a local (but 
uniform in space) mobility is invalid and detailed simu- 
lation of the transport of interacting carriers is unavoid- 
able. 

Effect of the interaction for the hopping transport in 
the DG model and GDM has been studied in recent paper 
[5^. It was found that the spatial correlation manifests 
itself even in the case of high carrier density. Indeed, a 
general tendency for the DG model is that transformation 
of the mobility curve with the increase of average fraction 
of the occupied sites p = na^ resembles the corresponding 
transformation of the curve with the increase of T (com- 
pare insets in Fig. [3 for example, with Fig. 1 in Ref. 



[27|): with the increase of p mobility becomes greater 
and the slope of n{E) curve becomes smaller. This is not 
the case for the GDM: here only mobility curve moves 
upward but the slope remains approximately constant. 

This difference could be easily understood. It was 
noted that the field dependence of ^ in the GDM is gov- 
erned by the carrier escape from deep states to the near- 
est sites having much higher energy, and the field-induced 
shift of site energies leads to Eq. (pn|) . Random charge 
distribution provides a smooth random energy landscape 
superimposed on the intrinsic disorder, but typical addi- 
tional variation of energy at the scale a is negligible for 
small p. Hence, estimation ([20)) remains valid and the 
slope of the mobility curve does not depend on p. 

Situation in the DG model is different: here mobility 
field dependence is governed by the carrier escape from 
critical clusters, as described in Section IIII Bl If we in- 
crease the density of carriers, then at first they fill these 
critical traps, because the release time is maximal here. 
Hence, transport of more mobile carriers is governed by 
clusters with the size that differs from (it is smaller). 
This means that the effective critical size rj;f depends 
on p. According to Eq. (|14l) . this is equivalent to the 
introduction of the effective temperature Toff, depending 
on p, and Tcff grows with p. This conclusion is in good 
agreement with Fig. [T] while for the GDM Teff does not 
depend on p and is very close to T, for the DG model 
Tcff monotonously grows with p. 

It was found also that for the DG model in the case of 
moderate p < 0.1 carrier drift mobility increases with p, 
exactly as in the case of non-interacting carriers ]^, [sl] . 
Yet modification provided by the interaction is still sig- 
nificant (see Fig. [8]). There is a striking difference be- 
tween the effects of interaction (i.e., carrier-carrier repul- 
sion) on the mobility in the DG model and GDM. In 
the DG model repulsion between carriers makes mobil- 
ity even greater than in the case of no interaction, while 
for the GDM the opposite situation takes place. This 
difference agrees with the cluster structure of the DG: 
if a carrier is trapped by some valley of the energetic 
landscape, then the whole valley with many sites having 
low energies becomes blocked for other carriers because 
of repulsion. Thus, filling of the deep states is much 
more effective in correlated landscape if carrier repulsion 
is taken into account. This is the reason for the increase 
of carrier mobility in DG, in comparison with the case 
of non-interacting carriers. No such effect takes place 
for the GDM, and here, evidently, the effect of charge- 
induced energetic disorder is responsible for the decrease 
of mobility in comparison to non-interacting case. 

This particular result disagrees with the result of a re- 
cent paper by Zhou et al [561 ]. They found that in the 
GDM carrier interaction enhances mobility in compari- 
son to the case of no interaction if a/kT ^ 1. This is op- 
posite to our findings. Quite probably, the disagreement 
stems from the under-relaxation of the initial (random) 
carrier configuration used in Ref. 56] ; the relaxation pro- 
cess is pretty slow for interacting carriers if a/kT 3> 1. 
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FIG. 8. Dependence of the carrier mobility on the effective 
strength of the charge-charge interaction Uc ~ /eaa for 
p = 0.008, kT/a = 0.3, and eaE/a = 0.1. Black squares show 
the data for the DG model and empty circles show results for 
the GDM. For a typical disordered organic material Uc — 5. 
Here /i(0) is the mobility for Uc = 0. 



We cannot make more detailed comparison because typi- 
cal relaxation times are not provided in Ref. [s^ (in fact, 
even the strength of carrier repulsion is not provided). 
Our data indicates, for example, that for a/kT = 4 re- 
laxation is not completely over even for t/tQ = 1 x 10^ 
(see Fig. [S]); at that time carrier has already traveled in 
the field direction the distance of ~ 4 x 10'^ a. 

Remarkable feature of Fig. |9] is a universality of the 
late relaxation stage. Very early relaxation is different 
for the initial random distribution and minimal energy 
distribution (where every carrier was placed at the site 
where the total energy, provided by the intrinsic disor- 
der and all previously added carriers, has a minimum), 
but after t/tQ ~ 10 relaxation curves merge into a single 
curve. 



B. Effects of interaction: comparison with 
experiment 

One can suggest that carrier transport in organic field- 
effect transistors (OFETs) should be a natural choice 
for comparison of the simulation results with experiment 
[60l [6l| . Estimation of the carrier density in OFETs 
show that the density as high as 3 x 10^^ cm^"^ could be 
achieved [62i|, that for a « 1 nm corresponds to p ~ 0.03. 
Experimental data for the particular OFET should be 
compared with the TOF data for a sandwich device hav- 
ing transport layer of the same material; in this way we 



-1.5 



-2.5 



-3.5 



0/kT = 4 
eaE/a = 1 
p = 0.0625 




-4 -2 2 4 6 8 10 12 
ln(t/to) 



FIG. 9. Relaxation of the mean carrier velocity in the GDM 
for a/kT — 4, eaE/a = 1, and p = 0.0625 (squares). Empty 
squares show relaxation for random initial locations, while 
filled squares show relaxation for the case, when initial loca- 
tions have been taken at the minimal energy positions. Cir- 
cles shoe the relaxation for higher temperature a/kT = 3.33. 
Time and velocity scales are to = l/vo and vo — a /to, corre- 
spondingly. 



could use transport characteristics (e.g., cr), relevant for 
the intrinsic disorder in the material. Quite frequently, 
OFETs demonstrate mobilities much higher that the mo- 
bilities measured in TOF experiments, and usually mo- 
bility increases with the increase of p [62] . This fact is in 
general agreement with the model studied in this paper. 

However, careful analysis reveals much more compli- 
cated situation. Indeed, in many aspects OFETs are 
very far away from the model, considered in the cur- 
rent study. First of all, in OFETs transport occurs in a 
thin layer, close to the gate insulator. Quite probably, 
especially in polymer devices, structure of this layer dif- 
fers from the structure of the same material in the bulk 
(polymer chains could be arranged in a special way at the 
gate insulator surface). This arrangement could provide 
more ordered structure with less degree of energetic dis- 
order, thus mobility should be enhanced near the inter- 
face, but accumulation of surface defects and impurities 
at the interface could lead to the decrease of the mobility. 
Next, there is a clear indication that the roughness of the 
organic semiconductor/dielectric interface affects carrier 
mobility [gI]. At last, the very nature of a gate dielec- 
tric (specifically, its polarity) affects carrier mobility in 
OFETs, because a random orientation of polar groups in 
the vicinity of a transport layer induces an additional en- 
ergetic disorder in semiconductor [63.[65|. In short, there 
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are a lot of reasons to believe that transport properties of 
OFETs are too complicated to be directly compared with 
the results of this study. We can only state that a sig- 
nificant increase of the carrier mobility with the increase 
of carrier density in carefully manufactured OFETs does 
not contradict the results of our study. 

VII. CONCLUSION: FUTURE DEVELOPMENT 

We may conclude that the spatially correlated elec- 
trostatic energetic disorder is a dominant factor directly 
dictating major features of the hopping charge transport 
and injection in amorphous organic materials. Strong 
spatial correlation is an inevitable property of the ran- 
dom energy landscape, created by long range electro- 
static sources. It is worth to note that this inevitability 
is still not recognized by a significant part of the commu- 
nity. Sometimes one can read in papers that "our analysis 
does not indicate that there is a need to assume a certain 
spatial energy correlation" (see, e.g. the paper [s^). Typi- 
cally, it is believed that the non-correlated distribution of 
random energies is something natural (probably, because 
of the simplicity of the conception), but the correlated 
distribution needs the specific reasons to occur. We hope 
that the discussion in this paper demonstrates that in 
organic glasses quite the opposite is true: they are corre- 
lated media by very nature. If, for some reason, there is 
a need to suggest a non-correlated energetic disorder for 
charge carriers, then the very existence of such disorder 
in organic glasses is very difficult (probably, impossible) 
to justify. 

At the same time, our knowledge of some important 
aspects of the hopping transport in correlated landscape 
is not sufficient. We may mention the almost absolute 
absence of reliable (all the more so, exact) analytical re- 
sults for the 3D transport, scarcity of theoretical results 
on the shape of the current transients, and other open 
problems. 

We would like especially emphasize the ultimate deficit 
of the experimental data on the local orientation order 
in organic glasses: the local order could change the cor- 
relation function for short distances and, hence, change 



the mobility field dependence for strong fields. Our cur- 
rent knowledge in this area can be estimated using re- 
cent papers |66l - l69l |. it is absolutely insufficient for the 
reliable consideration of transport problems. Study of 
the transport of interacting carriers is still in its infancy. 
Reliable experimental evidence for the modification of 
electrostatic disorder at the electrode was not provided. 

Several recent papers advanced a program of so-called 
multiscale modeling of charge transport in amorphous 
organic materials |70l - l72| . According to this program, 
the very structure of the organic glass is simulated, us- 
ing some variant of the molecular dynamics or related 
approach, then the relevant parameters (positions of en- 
ergy levels, hopping probabilities) are calculated for ev- 
ery particular particular realization of the structure of 
the glass, and then the Monte Carlo simulation of the 
hopping transport is carried out. Using this approach, 
it was possible in some cases to calculate the mobility 
value, pretty close to measured in experiments (the dif- 
ference is about one order of magnitude) [7i|. At the 
same time, typical calculated mobility field dependences 
are much less steep, than the experimental ones. Prob- 
ably, the reason for this difference is insufficient size of 
the basic simulation sample (about 10^ molecules). Obvi- 
ously, limitations of the multiscale simulation are mostly 
determined by the achieved computer performance which 
is constantly increasing. Hence, this approach seems to 
be a very promising line of future investigation. In addi- 
tion, it could provide a valuable information on the local 
structure of organic glasses. 

This short and, by no means, exhausting list of open 
problems and possible directions of further investigation 
could serve as a proof for the reader that the study of 
hopping charge transport in organic materials will be a 
vibrant area of research for the observable future. 
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